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The electronic structure of interfaces between LaAlOs and SrTiOs is studied using local spin 
density approximation (LSDA) with intra-atomic Coulomb repulsion (LSDA+U). We find that the 
nature of the interface metallic states is strongly affected by the type of the structure (sandwich or 
bilayer) and by the termination surface of LaAlOs. In all structures the atomic relaxation plays a 
crucial role in the electronic properties of the system. While in sandwiches the structural relaxation 
produces a significant polarization in SrTiOs and Jahn- Teller like splitting of Ti 3dxy orbitals, in 
A102-terminated bilayers the relaxation occurs primarily in LaAlOs and results in an insulator- 
metal transition which has been observed experimentally with increasing thickness of the LaAlOs 
layer. 
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I. INTRODUCTION 

Interfaces in complex transition metal oxide het- 
erostructures can play an important role for their phys- 
ical properties^. This is especially important for het- 
erostructures which contain nm-thick films deposited on 
structurally compatible substrate layers. The lattice mis- 
match and the charge polarity due to structural disconti- 
nuities at a film-substrate interface are the driving forces 
of interfacial reconstruction which changes dramatically 
the interface electronic state and possibly the physical 
behavior of the entire film. 

The character of such a reconstruction can range from 
purely electronic to mixed electronic-lattice. The elec- 
tronic reconstruction occurs through a compensation of 
the interface polarity by charge self-doping. In addi- 
tion, the lattice relaxation can modify the band struc- 
ture (Jahn- Teller like effect) and redistribute the charge 
between the electronic levels with consequent changes of 
the electronic properties of the system. 

Besides the interface region, the film surfaces fre- 
quently undergo a surface reconstruction. For the sur- 
faces of transition metal oxides, this reconstruction can 
be related to an electron depletion of (i-orbitals accom- 
panied by lattice contraction, a process observed on sev- 
eral surfaces like for instance at Pt(lll)^. The charge 
compensation at the surface can also proceed via a for- 
mation of oxygen vacancies which bring extra positive 
charge into the top layer. In both cases, the lattice re- 
laxation and polarization plays a key role and leads to 
changes of chemical bonding and covalencies of surface 
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An example of a film-substrate heterostructure, which 
exhibits a combination of surface and interface recon- 
struction, is a system of a LaAlOs film on top of a SrTiOs 
substrate. This bilayer attracted much attention recently 
due to the discovered transition between metallic and in- 
sulating properties^, possible interface magnetism^ and 
superconductivit}^. Electrostatic tuning of the charge 



carriers allows to switch between superconducting and 
insulating state^. A nanoscale control of the interfacial 
properties has been achieved^ and nanoelectronics has 
become feasible^. 

In LaAlOs-SrTiOs-heterostructures, a several unit-cell 
thick LaAlOs (LAO) film is grown on a mm-thick SrTiOs 
(STO) substrate. Due to the extremely small thickness of 
the LAO film, the top surface states are almost directly 
coupled to the interface. This coupling leads to a new 
cooperative mechanism of reconstruction which involves 
both surface and interface electron states. As a conse- 
quence of such a combination, one may expect a strong 
dependence of the physical properties on the thickness 
of the LAO film. In fact, experiments^ with LAO-STO 
give clear evidence of an abrupt increase of the carrier 
concentration and conductivity in heterostructures with 
LaO-Ti02 interfaces of n-type and AIO2 termination of 
the LAO film. These abrupt changes have a character of 
a transition between insulating and metallic states which 
occurs upon an increase of the LAO-thickness above a 
critical value of 3 unit cells in the direction perpendicu- 
lar to the interface. As the pure electronic compensation 
of the interface polarity would result in 0.5 electrons in 
the interfacial region and in metallic properties, the in- 
sulating character of the system with the ultrathin LAO 
film suggests that other reconstruction mechanisms, dif- 
ferent from the purely electronic, play a dominant role. 

In this work, we present the results of electronic struc- 
ture calculations of LAO-STO superlattices. The struc- 
ture of model superlattices contains the bilayers of LAO 
and STO with a slab of vacuum on top of the LAO layer, 
a structure which closely resembles the experimentally 
studied LAO-STO heterostructures. To understand bet- 
ter the surface reconstruction, we compare our results 
with calculations for a sandwich-type STO-LAO-STO su- 
perlattice where an LAO unit cell is directly connected to 
the Ti02-terminated STO layers. In the analysis of the 
lattice degrees of freedom, we perform a full structural 
relaxation of the considered systems which includes the 
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FIG. 1: Schematic view of relaxed LAO-STO sandwiches 
with A^LaO = 2 planes of LaO and A^Ti02 — 2 and 3 planes 
of Ti02 in each substructure. Here, A^Ti02 — 2 for the bot- 
tom and A^Ti02 — 3 for the top structure, respectively. The 
distinct buckling in the interface Ti02 planes is evident. 

optimization of the lattice constants and of the atomic 
positions in supercells. In our calculations, we employ 
the SIC variant of LSDA+f/— on grids containing from 
(8 X 8 X 1) to (9 X 9 X 1) /c-points where the additional 
corrections due to Coulomb repulsion on 3d orbitals of 
Ti {UTi = 2 eV) and bd orbitals of La (/7La = 8 eV) are 
included, according to estimates in Ref. |l3|. 

As appears from our analysis of the Ti02-LaO inter- 
face stack, the reconstruction in STO-LAO-STO sand- 
wiches differs from that in STO-LAO bilayers. The rea- 
son is the significance of the additional distortions in the 
near-surface planes of LAO in the bilayers, whereas in 
sandwiches we find only small atomic displacements in 
the interface LaO layers. Moreover, the reconstruction 
of STO-LAO bilayers depends on the type of surface ter- 
mination of the LAO film. Below we consider the conse- 
quences of the surface termination. We will see that, de- 
spite similar Ti02-LaO interface stacking, the structural 
relaxation in all these systems leads to qualitatively dif- 
ferent properties which demonstrates the key role of the 
lattice reconstruction. 



II. INTERFACE RECONSTRUCTION IN 
SANDWICHES 

It is well established in the literature^ that for the 
polar (LaO)+-(Ti02)^ interfaces, the "polar catastro- 
phe" can be prevented by an electronic reconstruction 
mechanism^^. In sandwiches LAO/STO/LAO/. . . with 
LAO terminated by a LaO plane (see Fig. [1]), electronic 
reconstruction results in 0.5 electrons per interface unit 
cell required to maintain the overall neutrality of the sys- 
tem. It is remarkable that the negative "extra charge" is 
located mostly in the Ti 3d orbitals ^i^ii^ 

In the analysis of heterostructures which contain 
SrTiOa with its incipient ferroelectric properties, one also 
has to take into account the lattice degrees of freedom. 
Due to the high polarizability of the titanates, the charge 



polarity at interface regions will inevitably induce local 
atomic distortions. The resulting atomic relaxations will 
modify the interface band structure and metallicity which 
can bring radical changes into the electronic reconstruc- 
tion mechanism. 

In the investigation of a combined electron-lattice re- 
construction in LAO/STO, we first optimized the in- 
terface distance between Ti02 and LaO layers within 
LDAii. The optimal distance (1.89 A) is shorter then the 
bond length [Ti-0] (1.953 A) in bulk STO which implies 
a contraction of the interfacial TiOe octahedra in agree- 
ment with the results of Ref. 17. After this optimization 
at the interface, the local positions of atoms in the sand- 
wich were relaxed to the new energetically favourable val- 
ues. In each plane, the structural relaxation leads to sig- 
nificant, diametrically opposed displacements of Ti and 
O atoms in the [001] direction (denoted below by z). As a 
consequence, a static polarization is built up, character- 
ized by the elongation of TiOe octahedra by about 0.1 A 
in [001] direction and by a buckling of the Ti02 planes. 
The changes of the length Aao = ^a — zq of each (A-0) 
bond (A = Ti, Sr, La) are shown in Fig. [2j These dis- 
placements cause a formation of static dipole moments in 
an initially paraelectric STO layer (similarly to what has 
been found in GGA and GGA+/7, Ref. 18) and reflect 
the degree of buckling of LaO, Ti02 and SrO planes. A 
good convergence of the dipole moments within the in- 
terfacial region of 25 A is achieved already for systems 
containing seven Ti02-planes in the supercell. The de- 
pendence of the displacements on the distance from the 
interface agrees well with the results of Ref. [l9|, which 
suggest an exponential decay of the displacements. The 
largest bond distortions of 0.15 A appear in the Ti02 
planes and increase for increasing thickness of the STO 
layers. In contrast to the STO layer, the distortions in 
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FIG. 2: Static dipole distortions in the planes SrO and Ti02 
which form with increasing number of STO unit cells in a pe- 
riodic STO-LAO-STO supercell. The local planar distortions 
(vertical axis) are the changes of the length Aao = za — zq of 
each (A-O) bond (A = Ti, Sr, La) with respect to their bulk 
values. Here the coordinate z measures the distance from the 
interface LaO plane of the sandwich LAO/STO/LAO to each 
plane in STO. 
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FIG. 3: Projected densities of states for Til in the inter- 
face Ti02 layer and for Ti3 in the Ti02 plane separated by 
a distance 2a = 7.81 A from the interface in a STO-LAO- 
STO sandwich with unrelaxed and optimized atomic positions 
(A^Tio^ =3and ArLaO = 2). 
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FIG. 5: Electron density profiles in structurally optimized 
STO-LAO-STO sandwiches with A^lao = 2 for different 

A^TiOs. 



the LAO planes are negligibly small. It is worthwhile 
to note that the most significant atomic displacements 
A > 0.05 A do not extend beyond a thickness of ap- 
proximately 20 A, which measures an effective depth of 
structural relaxations near the interface. 

In the electronic band structure calculated within 
LSDA+/7, the deformation of TiOe produces a Jahn- 
Teller splitting of t2g states. It should be noted that the 
t2g states are partially filled by the charge carriers ap- 
pearing due to the electronic reconstruction of the polar 
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FIG. 4: Electron density profiles in structurally unrelaxed 
and optimized STO-LAO-STO sandwiches with A^Ti02 — 3 
and A^LaO = 2. The density is determined from the integra- 
tion over the density of states of Fig. [3l 



LaO-STO interface. Specifically, the structural distor- 
tions split the t2g levels into d^y and degenerate dxz and 
dyz bands. As the atomic displacements strongly depend 
on the distance from the interface, the level splitting and 
related changes in the band structure will be different in 
the interface region and in the layers more distant from 
the interface. Fig. [3] shows the modifications of the pro- 
jected density of states of Ti in the layer next to the 
interface (Til) and in a more distant (3^^) layer (Ti3). 
In the interface layer, the relaxation leads to a slight shift 
of the conduction 3dxy band and, more prominently, to 
a reduced DOS close to the Fermi energy which suggests 
a decrease of the charge carrier density. In contrast to 
the interface Ti02, the relaxed 3dxy band of Ti3 is dis- 
placed further below the Fermi level and the DOS at the 
Fermi energy is enhanced. The shifts of the conduction 
bands also produce changes in their respective occupan- 
cies. They can be calculated through the integration of 
the DOS in the range between the bottom of the conduc- 
tion band and the Fermi level. The calculated profiles of 
the electron charge density (Fig. |4]) suggest a charge de- 
pletion of the interfacial Ti02-LaO stack and consequent 
accumulation of the electron charge in the Ti02 layers 
which are more distant from the interface. 

The depletion of conducting charge in the interface lay- 
ers not only results from the relaxation but may also be 
triggered through a growing number of Ti02 planes in 
STO. This fact can be observed in Fig. [5] which demon- 
strates a redistribution of a fixed amount (0.5) of com- 
pensated charge with an increasing number A^Ti02 
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FIG. 6: Electron charge as a function of the number of Ti02 
planes in structurally optimized STO-LAO-STO sandwiches 
with A^Lao 2. 



Ti02 layers from 2 to 5. An important consequence of 
the electronic reconstruction at the interface is a decrease 
of electron density in the Ti02 and LaO interface planes 
(see Fig. [6|) which is stronger in systems with wider STO 
substructure. In our calculations, an increase of the num- 
ber of Ti02 planes, A^Ti02 7 leads to a decrease of charge 
in the interfacial Ti02 plane, with a tendency for a sat- 
uration of electron charge in the interface Ti02 plane to 
a value of about 0.1 electrons for a growing number of 
Ti02 planes in the STO slab. In the STO layer with 
5 Ti02 planes, the residual 0.4 electrons are distributed 
within the more distant Ti02 planes, but about 73% of 
the whole charge is located in the first three Ti02 planes. 
It suggests a confinement of electron charge in an inter- 
face region of a thickness of 10-20 A in the STO layer. 
This finding agrees well with the experimental distribu- 
tion of the Ti excess charge discussed in Ref. 4, where 
the distinct peaks of the excess charge and its substantial 
decay with increasing distance from the interface demon- 
strate the confinement of charge to the interfacial region 
of a width of about 30 A. Specifically, for sandwiches 
containing 9 to 11 Ti02 layers, the Ti excess charge ap- 
proaches a value slightly above 0.1 electrons at a distance 
of about 3-4 A, an experimental value^ which is compat- 
ible with the values presented in Fig. [5l 

We reason that the small width of the interfacial re- 
gion is consistent with the conclusion that a quasi-two- 
dimensional electron gas is formed near the interface 
of a LAO-STO heterostructure with a mm-thick STO 
substrate (as discussed in Ref. 5 and, more recently, in 
Ref.i20,). 



III. ROLE OF SURFACE RECONSTRUCTION 
FOR THE ELECTRONIC PROPERTIES OF 
BILAYERS 

Experimentally, it has been established that LAO-STO 
heterostructures with more than three unit cells of LAO 



on STO substrates are metallic. For the analysis of the 
insulator-metal transition with increasing thickness of 
the LAO film, one has to account for the features inher- 
ent to surfaces of transition metals and their oxides. The 
current studies of these surfaces are motivated by drastic 
changes of their physical and chemical properties with 
the realization of different types of atomic reconstruc- 
tion in the surface and subsurface layers. Furthermore, 
a number of ab-initio studies of this type of interfaces 
have identified a metallic n-type state with charge carri- 
ers which predominantly fill 3dxy orbitals.^^^^^ 

Recent experiments give indications of a temperature 
dependence of the surface termination of LAO which 
changes from AIO2 to LaO planes with an increase of 
temperature above 523 Ki^i^i^ Therefore, to under- 
stand the interface reconstruction mechanism in ultra- 
thin LAO films with an interference between surface and 
interface states, it is instructive to consider separately 
two possible types of the surface plane (AIO2 and LaO) 
in STO-LAO bilayers. 

For our theoretical studies, we have generated super- 
cells which contain several LAO unit cells deposited on 
1.5 SrTiOs unit cells. We have restricted calculations 
to structures with A/'Ti02 — ^ due to clear indications 
of charge confinement to a few unit cells of STO found 
in the analysis of sandwiches (section (III). To consider 
the surface properties, we have also introduced a 13 A 
vacuum layer on top of LaAlOa. The thickness of the 
LaAlOs film has been varied from 1 to 3 unit cells. In a 
system with the topmost plane LaO, an extra half of a 
LAO unit cell has been attached to the top of the LAO 
film. 

The LSDA+/7 calculations result in a stable nonmag- 
netic state of the system. In this state, the electronic 
properties strongly depend on the structural relaxation 
which is controlled by the thickness of the LAO film 
in the slab. Specifically, in the LSDA+/7 calculations, 
the initially chosen configurations with nonzero magnetic 
moments and charge inhomogeneities rapidly decayed to 
zero magnetic moments and to equal charge for all Ti 
atoms in each Ti02 plane which excludes the possibil- 
ity for a decisive role of charge or magnetic order on the 
resultant electronic state. 



A. AIO2 termination 

To study the bilayers STO-LAO with an AIO2 plane 
at the top surface (Fig. [7]), we have generated a num- 
ber of supercells with A^LaO varying from 1 to 3 LAO 
unit cells. In these supercells, the structural optimiza- 
tion results in a structure with large opposite distortions 
of cations and anions in the LaO planes and only small 
atomic displacements in the STO substrate (see Fig. [7]). 
Such a relaxation pattern is completely distinct from the 
systems with LaO termination which will be considered 
in the following section IIIIBI The details of the struc- 
tural changes are presented in Table HI 



TABLE I: Atomic displacements (in A) in the top planes 
of STO-LAO bilayers with an AIO2 plane at the surface. 
Here the atomic displacements of the top La sites are defined 
as A^La = ^La — ^La whcrc z^a IS the planar unrelaxed 
coordinate, and of the oxygen sites as Azq = zq — Zq. 
An additional index La or Al at 2;oLa/Ai refers to the 
oxygen position in the La or Al planes, respectively. The 
atomic distortions in the top LaO plane are defined as 
AzLaO = BAzLa — 2Azoj^^ The buckling of the surface is 
given by the atomic distortions in the top AIO2 layer which 
are defined as A^aio2 — ^A^Oai ~ 3Azai- 



A^LaO 


A^La 




AzLaO 


AZA\ 


Azoai 


A2;aio2 


1 


0.29 


0.05 


0.77 


0.015 


-0.02 


-0.125 


2 


0.11 


-0.07 


0.46 


-0.11 


-0.13 


-0.19 


3 


0.06 


-0.09 


0.36 


-0.09 


-0.09 


-0.09 



The first important feature is the strong buckling of 
the LaO planes where the relative [La-0] displacement 
A^LaO = 3zLa " 2zo approaches a maximal value of 
0.77 A in the structure with A^Lao = 1 (cf. Fig^. This 
displacement parameter A^LaO is introduced so as to 
characterize the dipolar distortion of a LaO plane. For 
A^LaO > 1 the displacement parameter AzLaO is consid- 
erably reduced. 

The second key property is a displacement of the top- 
most AIO2 plane. It comprises significant inward dis- 
placements of both Al and oxygen towards the LaO plane 
and a consequent decrease of the distance between the 




FIG. 7: Schematic view of relaxed LAO-STO heterostruc- 
tures with A^Ti02 — 2 and with top layers of LaAlOs ter- 
minated by AIO2. Here the configurations with A^LaO = 1 
and A^LaO = 2 are shown. For A^LaO = 1, the covalencies 
between La and the titanium and oxygen ions in the plane 
below are weak and no bonds are drawn. The full supercell 
for our calculations contains, apart from the vacuum cells, 
also the inverted structure so that the supercell geometry has 
two inversion symmetric interfaces and surfaces. 
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FIG. 8: Electron densities of states of a relaxed STO-LAO 
bilayer with an AIO2 plane as the surface plane for A^xiOs — 2 
and different A^LaO- 



AIO2 plane at the surface and its adjacent LaO plane. 

In the case of A^LaO = 1, the electric dipole moment, 
which is induced by the atomic distortions, is given by 
APlao = Qi^Zi where Qi refers to the charge of 
the i-th site {i = Al, La, Oai, OLa) where OAi/La is 
the oxygen in the A102/LaO plane. One verifies that 
APlao = ^lAzLaO + Z2AZA102 = 0.9 eA where Zi = 
—Z2 = e. Here, the displacement parameter for the AIO2 
plane is AzaiOs = 4A2:oai - SAzai- 

On the other hand, the polar catastrophe in 
LAO^^ leads to the net electric polarization Plao — 
— eaLAo/2 = —1.9 eA, where aLAO — 3.8 A is the ex- 
perimental bulk lattice constant of LaAlOa. While the 
interaction of electrons with the polarization Plao t^^^ds 
to close the energy gap and stabilize the metallic state, 
the polarization due to atomic distortions is of opposite 
sign with respect to Plao leads to a competing ten- 
dency to stabilize the insulating state. In the monolayer 
LAO film with a dielectric constant £, the interaction 
energy E = ePLAo/(47r6:o^<^LAo) ^ — eV of the to- 
tal dipole moment Plao = ^lao + APlao = — 1 eA 
with the electronic charge results in a reduced but finite 
energy gap of 1.5 eV which, for comparison, is approxi- 
mately 4 eV in the LDA calculations of bulk LAO. 

For the structure with A^LaO = 2, we have ^ao ~ 
— ^LAoeaLAo/2 = —3.8 eA and APlao = 0.65 eA. As a 
result, the increased net polarization Plao = —3.15 eA 
closes the energy gap and leads to the metallic state ob- 
served in Fig. [HI 

In the densities of states in Fig. [HI the buildup of metal- 
lic states proceeds through the appearance of the electron 
carriers in the Ti 3dxy bands which is exhibited in Fig.[TQl 
This is accompanied by the filling of the gap which, for 
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^LaO = 1, exists between the oxygen 2p and the La 5d 
and Ti 3d bands. The 2p states of the oxygens in the 
AIO2 planes approach the Fermi energy from below and 
the interface Ti 3dxy states cross the Fermi energy from 
above (see Fig. [9]). 

In experiments, STO-LAO bilayers with LAO films of 
^LaO = 3 are still insulating whereas LAO films with 
A/^LaO ^ 4 have been found to be metallic. In the present 
LDA-calculations, due to the underestimation of the en- 
ergy gap by LDA, the "critical thickness" of the LAO 
layer is decreased by 2 and the metallic state occurs al- 
ready at A^LaO = 2. Moreover, other localization effects 
may play a role in real materials^^^^ , especially when 
the density of states at the Fermi energy is still small (cf. 
Fig. [9]). In GGA calculations without lattice relaxation, 
the transition to a metallic state was identified for LAO 
films with A^LaO = 4 (supplement of Ref. 9). The GGA 
evaluation with lattice relaxation^^ produces a transition 
at around five monolayers of LAO. 

Below the critical LAO-thickness, the polarity compen- 
sation via self charging (see Chap. ^ reduces the energy 
gap, a process which however is partially suppressed due 
to structural distortions in the LAO layer. Similarly to 
the GGA results presented in Ref. 113, we can conclude 
here that the polarization in the LAO layer due to struc- 
tural relaxation competes with the electric field caused 
by the polarity catastrophe and in this way establishes 
the insulating state in STO-LAO bilayers below the LAO 
critical thickness. In contrast, above the critical thick- 
ness the polarity compensation through the surface and 
interface charging leads to the metallic properties. 

In the considered bilayers, two different charging pro- 
cesses can be distinguished which are exhibited by the 
composition of the density of states close to the Fermi 
energy (see the two panels of Fig.[9j). Near the positively 
charged LaO layer at the Ti02-LaO interface, the com- 
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FIG. 9: Electron densities of states of Ti and oxygens in 
the LAO layer of a relaxed STO-LAO metallic bilayer with 
A^TiOo 2 and A^LaO 2. 
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FIG. 10: Projected electron densities of states of the interface 
Ti 3dxy orbitals for relaxed STO-LAO metallic bilayers with 
^Ti02 — 2 and A^LaO = 1, 2 and 3 from bottom to top. 



pensation produces the additional negative charge accu- 
mulated at Ti 3dxy orbitals (see Fig. [9] and Fig. [TQ|). In 
distinction to the interface, the negative polarity near the 
surface AIO2 plane leads to additional holes distributed 
over the 2p orbitals of the surface oxygen which is clearly 
observed in the bottom part of Fig. [9l Recently, indi- 
rect indications for positive surface charge in STO-LAO 
have been discussed in Ref. 9 where such a positive sur- 
face charge is compensated by oxygen vacancies. In other 
experimental conditions, the surface hole charge can be 
also neutralized by surface oxidation or by near-surface 
bonding with negatively charged radicals which can be a 
possible reason for the insulating state usually found on 
the AIO2 surfaces of LAO-STO bilayers. 



B. LaO termination 

The low- and room-temperature experiments are per- 
formed with heterostructures, in which the LAO films 
are terminated by AIO2 planes in the case of Ti02-LaO 
interfaces. There are, however, several experimental in- 
dications of a transformation of the termination plane 
from AIO2 to LaO, which occurs at temperatures around 
423 K— Such a transformation is explained through 
a formation of oxygen vacancies with simultaneous dis- 
placements of Al towards the subsurface layers. Although 
the discussed AIO2-LAO transformation remains still an 
open question, it is highly instructive to consider LaO 
terminated LAO films in the bilayers STO-LAO and to 
compare the electronic states with the electronic proper- 
ties of AIO2 terminated layers. 



TABLE II: Structural parameters (in A) for different A^LaO 
in STO-LAO bilayers with LaO termination. Here the 
atomic displacements of the top La sites are defined as 
A;2;La = zi^a — ^La whcrc z^a IS the planar unrelaxed coor- 
dinate, and of the oxygen sites as Azq = zq — Zq. The 
buckling of the surface is given by the atomic distortions of 
the top LaO bond A^LaO = ^zi^a — '^zq, and the length of 
the LaO-bond is denoted by ALaO- The distance between 
the top LaO and the nearest subsurface plane is S. 



A^LaO 


a 


c 


6 


AzLa 


Azo 


A^LaO 


ALaO 


2 


3.805 


32.17 


1.8 


-0.12 


0.02 


-0.4 


2.7 


3 


3.755 


39.75 


L85 


-0.05 


0.05 


-0.25 


2.69 


4 


3.759 


47.32 


L85 


-0.05 


0.05 


-0.25 


2.66 



For LaO termination at the topmost surface of the bi- 
layers, we suggest and elaborate below that the recon- 
struction is controlled by the surface tensile stress which 
causes a charge occupation of 5d orbitals and a simulta- 
neous contraction of the top surface layers^. 

In our supercells shown in Fig. [TTl the thickness of 
the LaAlOs film has been varied from 1.5 to 3.5 unit 
cells whereas the thickness of the SrTiOa slab was fixed 
to 2 Ti02 planes. In the heterostructure, the po- 
lar (La0)+^/Ti02 interface results in an effective elec- 
tron doping by exactly one electron which is needed 
to achieve charge neutrality. In the context of the re- 
cent experiments^'^ the key question is the influence of 
the structural relaxation and the concomitant electronic 
reconstruction on the interface state, especially on its 
metallic and magnetic properties. 

For the optimization of such supercells, we have per- 
formed a relaxation of the lattice constants a = b in the 




FIG. 11: Schematic view of relaxed LAO-STO heterostruc- 
tures with A^TiOs — 2 and with 1.5 and 2.5 top unit cells of 
LaAlOa (A^LaO = 2 and 3, respectively). 
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FIG. 12: Electron density map in the {x,y) plane of (a) 
unrelaxed and (b) optimized STO-LAO bilayer with A^LaO = 2 
and A^Ti02 — 2. The coordinate z/c = 0.29 corresponds to 
the top surface LaO plane in the LAO film. 



(x^y) plane and of the subsurface interplanar distances 
jointly with a further relaxation of the local atomic posi- 
tions in the heterostructure. The results are presented in 
Table [ni As appears in our calculations, the optimized 
planar lattice constant, a, strongly depends on the thick- 
ness of the LAO film, varying from 3.805 A for Ni^^o = 2 
to 3.76 A in a supercell with A^LaO = 4 which clearly 
shows a tendency to approach the bulk lattice constant 
for LaAlOs. It is worth noting that the volume LDA- 
optimization of bulk LaAlOs gives ^lao ~ 3. 75 A which 
underestimates the experimental value aLAO = 3.789 A. 
From Table [III finds that already in ultrathin LAO 
films with A/'LaO = 3 and A^LaO = 4 the planar lattice 
parameters a are sufficiently close to the bulk values. 

In contrast to the sandwich-type structures character- 
ized by small atomic displacements in the LAO film, in 
the considered bilayer the atomic relaxation leads (i) to 
a substantial change of the distances S between the top- 
most LaO and nearest subsurface AIO2 layers and (ii) to 
significant opposite displacements of the upper La and 
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X (Angstrom) 

FIG. 13: Electron density map in the {x,y) plane of opti- 
mized LAO-STO bilayer with iVLaO = 4 and A^TiOs = 2. The 
coordinate z/c = 0.34 corresponds to the top surface LaO 
plane in the LAO film. 



oxygen atoms of the top LaO surface layers. Specifically, 
we find that the distance [LaO-A102] decreases by a value 
which ranges from 0.05 A in a system with A/'LaO = 4 to 
0.1 A for a structure with A/'LaO = 2 (see Table HI]) . Fur- 
thermore, for larger A^LaO the decrease of the La-0 dipole 
moments in z direction and of the contraction of the pla- 
nar lattice constant results in a shortening of length of 
the planar [La-0] bond which varies from ALaO = 2.7 A 
for A^LaO = 2 to ALaO = 2.66 A for A^LaO = 4. We also 
note that in an unrelaxed heterostructure, where the lat- 
tice constant a is fixed to asxo = 3.905 A, the length of 
the [La-0] bond is substantially longer than that of the 
relaxed system. To apprehend the role of the lattice and 
bond contraction for the electronic properties of the sub- 
surface layers, we performed LDA+/7 calculations of the 
electronic structure of the generated optimized supercells 
for different values of A^LaO- 




Number of LAO unit cells 

FIG. 15: Centers of gravity of top La 5(1^2 _y2 and interface 
Ti 3dxy orbitals versus A^LaO. 



As compared to the unrelaxed system, in a structure 
with a 1.5 unit-cell LAO film (A^LaO = 2) on SrTiOa, the 
Jahn- Teller like splitting due to atomic distortions leads 
to an increase of the gap between O 2p bands formed by 
oxygens in Ti02 and LaO planes and mixed La 5d and 
Ti 3d states, an effect similar to that in the STO-LAO- 
STO-sandwiches. An opposite trend is obtained when 
we increase the thickness of the LAO film. Specifically, 
the thicker LAO layers exhibit smaller energy gaps and 
consequently higher charge carrier densities in the top 
LaO layer. The structural relaxation plays a central role 
for the modifications of the electronic structure. This 
is also illustrated in the charge density maps shown in 
Fig.inja) and Fig.IHJb). 

In the unrelaxed structure with A^LaO = 2, more charge 
(about 0.6) is localized on La (Fig.[T2fa)) of the topmost 
LaO layer. In the optimized system, the lattice contrac- 
tion leads to a strengthening of covalent bonds and to 
a substantial carrier density which is distributed in the 
space between La and oxygen ions (Fig. [I2fb)). In the 
optimized systems with A^LaO > 3, the increasing con- 
traction results in a further increase of the charge density 
which is redistributed between La and O in the top LaO 
layer (Fig. [T3|) . The obtained feature is closely related 
to the changes of the electron occupancy of top La and 
interfacial Ti shown in Fig. [T4j As compared to the film 
with A/'LaO = 1, in a thicker LAO film with A^LaO = 3 
we find a transfer of the carrier density from the inter- 




Energy (eV) 

FIG. 14: Interface electron charge on La in the top LaO 

and on Ti in the interfacial Ti02 planes for different numbers FIG. 16: Partial density of states of La in the top LaO plane 
A^LaO of LAO planes in the heterostructure. for different A^LaO- 
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facial Ti02 planes to the top LaO layers which appears 
due to the described structural transformations in sub- 
surface layers. In Fig. [15] this Jahn Teller effect can be 
clearly detected by the strong decrease of the center of 
gravity of the 5dx2_y2 band of the top La which con- 
trasts with the interfacial Ti 3dxy band center that re- 
mains almost unchanged with the variation of the LAO 
thickness. Here a gravity center Eg of a band with a 
density of states p{E) is calculated from the expression 
Eg = J dEEp{E)/ J dEp{E). The consequential increase 
of the electron charge on La (Fig. [T4|) is also related to 
the shift of the bottom of the La bd^^-y^ band towards 
lower energies in Fig. [16] and is the direct result of the 
discussed electronic reconstruction. 



IV. CONCLUSIONS 

We have studied the electronic structure of interfaces 
between the bulk insulators LaAlOs and SrTiOa using 
the local spin density approximation (LSDA) with intra- 
atomic Coulomb repulsion (LSDA+/7). We find that the 
nature of the interface metallic states is strongly affected 
not only by the type (sandwich or bilayer) of the struc- 
ture but also prominently by the termination surface of 
LaAlOa. In all types the atomic relaxation plays a crucial 
role for the electronic properties of the system. 

In STO-LAO-STO sandwiches the structural relax- 
ation results in an elongation of the TiOe octahedra by 
about 0.1 A in the [001] direction and in a buckhng of the 
Ti02 planes of up to 0.15 A for the plane nearest to the 
STO-LAO interface. These displacements cause a signif- 
icant polarization in SrTiOa within about 10 A from the 
interface. They also result in a Jahn- Teller splitting of 
the i2g states and the interface charge resides mostly in 
the Ti ?tdxy bands. The conduction charge is confined 
essentially to the first three Ti02 planes. 

In bilayers the surface presents an additional impor- 
tant actor which may control the electronic state at the 
interface, especially when the surface is separated by an 
ultrathin film from the interface. For a LAO film on top 
of a STO substrate, the termination (either AIO2 or LaO) 



at the surface is crucial for the atomic relaxation and the 
electronic state at the LAO-STO interface. 

In A102-terminated bilayers the relaxation occurs pri- 
marily in the LaAlOa film with minor displacements in 
SrTiOa. Specifically, the buckling of the subsurface LaO 
plane and the modified distances between the topmost 
(oppositely charged) AIO2 and LaO planes introduces a 
polarization which decreases with the number A/LaO of 
LAO unit cells. For the unit cell thick LAO films, it is 
this lattice relaxation with the corresponding polariza- 
tion of the film which counteracts the polarization in- 
troduced through the self charging (by 1/2 electron to 
compensate the interface polarity). For films with one 
unit cell LAO we identify an insulating state character- 
ized by an enegy gap of 1.5 eV between the Ti Zdxy states 
of the Ti02 plane next to the interface and the oxygen 
states of the topmost AIO2 plane. For films with more 
than one unit cell LAO (i.e., A/'LaO > 1) this gap closes 
due to the enhanced polarization from the self charging. 
Then the metallic state at the interface is formed by the 
Ti ?>dxy band with a finite density of states at the Fermi 
energy. This scheme results in a transition from insulat- 
ing to metallic states: in the present LDA-calculations, 
due to the fact that LDA underestimates the gap, the 
metallic state occurs already at A^LaO = 2 whereas ex- 
periments show that bilayers with A/'LaO = 3 are still 
insulating whereas LAO films with A^LaO ^ 4 have been 
found to be metallic. 

Finally, for LaO-terminated bilayers the lattice re- 
laxation is different: the surface tensile stress causes a 
charge occupation of La bd orbitals and a corresponding 
contraction of the top surface layers. Here, one electron 
per interface unit cell is required to maintain the overall 
neutrality of the system. Both, the Ti '^dxy states in the 
plane next to the interface and the La bdx^ -y^Yidj^e finite 
density of states at the Fermi energy and stay metallic, 
irrespective of the number unit cells in the LAO film. 
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